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Abstract. —A minimum viable population (MVP) size is estimated for the critically endangered black-footed ferret 
by examining five basic methods: experiments, biogeographic patterns, theoretical models, simulation models, and 
genetic considerations. Each method is evaluated for its applicability to the ferret and endangered species in general 
with two criteria in mind: (1) potential research impacts to target species and (2) the value of scientific accuracy and 
precision in relation to short-term conservation needs. For the black-footed ferret, the genetic method proved to be the 
most useful, resulting in an MVP estimate of about 200 ferrets for maintenance of short-term fitness. For most 
endangered species, a combination of the simulation and genetic methods will probably yield the best estimate of MVP 
size. MVP estimates have direct implications for future research, management, and recovery efforts of endangered 
species. 


Following the 1981 discovery of the critically 
endangered black-footed ferret (Mustela ni- 
gripes) in northwestern Wyoming, the only pop¬ 
ulation of the species known, two primary goals 
of a ferret study were established: (1) conserva¬ 
tion of the population and (2) recovery of the 
species (Clark 1984a). Estimation of population 
parameters such as distribution, birth rate, 
death rate, and immigration/emigration was a 
major objective. In addition to providing valu¬ 
able life history data for the ferret, estimating 
these parameters would enable us to address a 
key question: What is the minimum ferret popu¬ 
lation size necessary for the species to persist 
(i.e., avoid extinction)? Because of the apparent 
small size of the Wyoming population in 1981, 
the answer was paramount for successful conser¬ 
vation and, ultimately, species recovery. 

The minimum population size concept for spe¬ 
cies conservation was recently introduced by 
Shaffer (1981), with potentially significant impli¬ 
cations for. endangered species programs. Al¬ 
though many endangered species studies have 
been conducted, few have tried to estimate min¬ 
imum population sizes (Shaffer 1978, 1981, La- 
cava and Hughes 1984, Salwasser et al. 1984). 
Although US DA Poorest Service regulations re¬ 
quire that minimum viable populations of en¬ 
dangered species be maintained on Forest Ser¬ 
vice lands, these population numbers are not 
estimated but are instead routinely established 
as those levels specified in recovery plans 


(Lacava and Hughes 1984). In this paper, we 
estimate minimum ferret population numbers 
and, as a corollary, minimum area require¬ 
ments. Our purpose is to establish guidelines for 
future ferret research and management and to 
provide impetus, direction, and encouragement 
to other endangered species programs to use the 
concept of minimum population size. 

Conceit of Minimum Viable 
Population Size 

The notion of a required minimum population 
size for species conservation was first embodied 
by Shaffer (1981) in his concept of minimum 
viable population (MVP) size. He defined the 
MVP for a species as the smallest isolated popu¬ 
lation having a 99% chance of remaining extant 
for 1000 years despite foreseeable effects of four 
types of stochastic events: (1) demographic 
stochasticitv, or the chance events in the sur¬ 
vival and reproductive fitness of a finite number 
of individuals; (2) environmental stoehasticity, 
or perturbations due to habitat parameters, 
competitors, predators, and disease; (3) natural 
catastrophes, such as randomly occurring floods, 
fires, droughts, etc.; and (4) genetic perturba¬ 
tions resulting from changes in gene frequency. 
Shaffer (1981) stressed the tentative nature of his 
definition and emphasized the importance of 
defining MVP with explicit but flexible criteria 
such as time frame and survival probability. For 
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example, the survival probability could be set at 
95% or any other level, and the time frame could 
be shortened or lengthened as appropriate for 
planning needs. 

Shaffer (1981) proposed five basic methods for 
determining MVP size: experiments, biogeo¬ 
graphic patterns, theoretical models, simulation 
models, and genetic considerations. In the next 
section, we elaborate on the application and util¬ 
ity of each method for estimating MVP size for 
ferrets in particular and endangered species in 
general. The usefulness of each method is prag¬ 
matically evaluated relative to: (1) acceptable 
levels of potential research impacts on target 
species and (2) the real value of scientific accu¬ 
racy and precision in relation to short-term con¬ 
servation needs. In endangered species recov¬ 
ery work, these two concerns must be addressed. 

The MVP concept is a rapidly evolving one. In 
late 1984, the Forest Service and other agencies 
sponsored a workshop to provide a state-of-the- 
art review of the MVP issue (M. L. Shaffer, 
personal communication). A major purpose of 
the workshop was to consider revisions of Forest 
Service procedures for estimating MVP size. 
The primary focus of the workshop was a discus¬ 
sion of the relative importance of demographic, 
environmental, and genetic variability in deter¬ 
mining minimum population numbers as well as 
the development of new simulation and analyti¬ 
cal models for estimating MVP size. Proceedings 
of this workshop, which will be published in 
1986, may provide additional methods not dis¬ 
cussed in this paper and/or may suggest revi¬ 
sions of methods outlined and discussed below 
for estimating MVP sizes of endangered species. 

Determining MVP Size 
Experiments 

This method establishes isolated populations 
of species and monitors their population dynam¬ 
ics through time. If censusing can be accom¬ 
plished without stressful manipulation of indi¬ 
viduals (as can occur in a capture/recapture 
program), this approach requires only low levels 
of research activities. The only necessary data 
are population numbers and the area inhabited 
over a specific time period (i.e., time frame in 
MVP definition). 

Determining MVP via experimental methods 
brings up three concerns. Although methodo¬ 


logically simple, this approach is not plausible 
for many species, particularly endangered ones. 
First, for many endangered species, there exists 
no possibility of studying several isolated popu¬ 
lations. In the case of the ferret, there is only one 
known population, and consequently there is no 
way to measure variability in persistence of dif¬ 
ferent-sized, isolated populations. Second, time 
is not an unlimited resource in most endangered 
species studies. The years necessary to monitor 
persistence of populations for MVP estimates 
are simply not available for most endangered 
species. Lastly, many wildlife species cannot be 
censused sufficiently without some type of mark/ 
recapture program. These techniques can have 
negative research impacts (e.g., trap mortality) 
that some endangered species projects may not 
be willing to risk, particularly in the early stages. 
Because endangered species programs in gen¬ 
eral possess high levels of uncertainty, there 
may be hesitancy to increase uncertainty to even 
greater levels (Clark 1984b). Such has been the 
case in the Wyoming ferret project, where one 
research team chose initially not to employ 
mark/recapture and radio telemetry procedures 
until other methods proved the population large 
enough to withstand some stress (Clark 1981). 

The experimental method could prove useful 
when extensive population data are available 
prior to a species becoming endangered. For 
example, mountain goat (Oreamnos ameri- 
canus ) and bighorn sheep (Ovis canadensis) 
populations have been reintroduced to several 
sites in the western United States in recent years 
(Rideout 1978, Wishart 1978). If wildlife re¬ 
searchers monitor the dynamics of these intro¬ 
duced populations, the data should eventually 
allow them to estimate MVP sizes for these spe¬ 
cies. Such information would be useful not only 
for future translocations of these species’ popula¬ 
tions but also for management if, for example, an 
isolated population showed a decline in num¬ 
bers. 

Biogeographic Patterns 

By studying the distribution patterns of spe¬ 
cies occupying patchy or insular habitats, an esti¬ 
mate of MVP size and minimum area require¬ 
ments can be obtained. However, populations 
examined should be in equilibrium and their 
approximate period of isolation known (Shaffer 
1981). Under these conditions, researchers can 
estimate the smallest area inhabited by a species 
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Fig. 1. Distribution of 37 white-tailed prairie dog colonies on the Wyoming black-footed ferret study area, 
1981-1984. 


and the percentage of patches of a particular size 
that a species occupies. The smallest population 
that has persisted over an extensive time period 
can then provide a first estimate of MVP size. 

Because the ferret occupies a patchy habitat 
(i.e., prairie dog colonies with relatively discrete 
boundaries. Fig. 1), this approach can be ap¬ 
plied, at least in part, to the Wyoming popula¬ 
tion. Old time trappers have reported catching 
ferrets in the 1920s and 1930s on the Wyoming 
study area, thus indicating that ferrets may have 
occupied the site for at least 50 years. Extinction 
and recolonization, although unlikely, could 
have occurred during this period. Population 
censuses over the last four years (Clark 1986) 
suggest that the Wyoming ferret population is 
not in equilibrium but is increasing. As a result, 
more sophisticated census techniques (e.g., 
mark/recapture) are being used, leading to more 
precise and accurate population data. In the fu¬ 
ture, researchers should thus be able to deter¬ 
mine with some degree of confidence whether 
the Wyoming population is reaching an equi¬ 
librium density. 

If the Wyoming ferret population stabilizes in 
numbers, then data from this population will 
provide the first field assessment of MVP size 


and minimum area requirements for the species. 
However, due to the time required to obtain 
these data, the biogeographical method will not 
likely be useful for short-term management de¬ 
cisions concerning MVP size. Nevertheless, 
useful information on minimum area require¬ 
ments can be gleaned from the biogeographical 
approach. 

In the short run, researchers will need to de¬ 
cide what size of an area is necessary for ferret 
reintroduction or translocation for recovery 
planning (Forrest et al. 1985; Houston et al. 
1986). A recommendation for an area larger than 
the Wyoming site is warranted for two reasons. 
First, two habitat patches (e.g., Wyoming and 
South Dakota) of equal size may not be the same 
in habitat quality. For example, Wyoming fer¬ 
rets occur on white-tailed prairie dog (Cynomys 
leucurus) colonies, w'hereas a South Dakota 
population of ferrets occupied black-tailed 
prairie dog (C. ludovicianus) colonies (Hillman 
and Clark 1980). White-tails usually form small, 
sparsely populated colonies; black-tails form 
large, densely populated colonies (Hoogland 
1981). Consequently, a habitat patch that is large 
enough in one part of the range of the species 
might be insufficient in another due to differ- 
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ences in habitat quality. Second, nothing is 
known about the frequency with which ferret 
populations go extinct on habitat patches of vari¬ 
ous sizes. Therefore, it would be inadvisable to 
rely on the size of the Wyoming study area as an 
exemplary model of area requirements for black¬ 
footed ferrets. 

Although the biogeographical approach can¬ 
not be used for endangered species with con¬ 
tiguous distributions, it is applicable for species 
with insular or patchy distributions. Its greatest 
use should be for endangered species that oc¬ 
cupy patches of different size and similar quality 
and for which data on population numbers and 
isolation periods are available. For example, the 
Shoshone sculpin ( Cottus greenei) is restricted 
to different-sized spring systems along a 45-km 
stretch of the Snake River in south central Idaho 
(Wallace et al. 1984). Relative and absolute den¬ 
sities of the sculpin in many of these spring sys¬ 
tems were determined during a status survey in 
1980-1981 0- Griffith, personal communica¬ 
tion). By monitoring the persistence of these 
spring populations of this short-lived species 
over the next several years, researchers could 
estimate MVP size and minimum area require¬ 
ments for this species. 

Theoretical Models 

Although there are several theoretical models 
that predict extinction probabilities and times 
for a population, most are too complex for the 
simple data bases of endangered species. The 
theory of island biogeography (MacArthur and 
Wilson 1967), however, has received consider¬ 
able attention in the conservation field. Though 
originally developed in relation to true oceanic 
islands, this theory has been widely applied to 
continental habitat “islands” or patches (e.g., 
Brown 1971). The theory’s greatest potential ap¬ 
plication has been in the design of nature pre¬ 
serves for maintenance of species diversity (for 
review, see Margules et al. 1982). Much less 
attention has been devoted to those aspects of 
the theory that predict the distribution of a sin¬ 
gle, insular species (Smith 1974, Fritz 1979). 

Because the Wyoming ferret population is lo¬ 
cated on a large habitat “island” (i.e., prairie dog 
colony complex in Fig. 1), we employed the T k 
model of island biogeography to predict proba¬ 
bilities of successful colonization and times to 
extinction for the ferret. Regarding the Wyo¬ 
ming prairie dog colony complex as a single is¬ 


land is supported by three lines of evidence: (1) 
several prairie dog colonies at the Wyoming site 
contained at one time or another only one adult 
ferret, suggesting that interbreeding of ferrets 
among prairie dog colonies is occurring, (2) 
movement data reveal that ferrets range an aver¬ 
age of 2.5 km between colonies, with the maxi¬ 
mum intercolony movement being 5.7 km (For¬ 
rest et al. 1985), and (3) aerial and ground 
surveys indicate that potential ferret habitat de¬ 
clines significantly beyond the boundaries of this 
single colony complex, strongly suggesting that 
these ferrets represent a distinct, isolated popu¬ 
lation. 

Data required in the T k model are carrying 
capacity (K), per capita birth rates (A), and per 
capita death rates (|jl ). Direct observations of 
ferret litters in Wyoming indicated an average 
litter size of 3.4, or 1.7 female young per adult 
female (Forrest et al., in manuscript). We as¬ 
sumed that the sex ratio approximated unity, 
natality did not vary with age, and all females 
bred. No data on juvenile mortality were avail¬ 
able for ferrets; data from other mustelids indi¬ 
cated a medium (60%) to high (80%) rate of juve¬ 
nile mortality (King 1980, M. nivalis , King 1983, 
M. erminea ; K. C. Walton personal communica¬ 
tion, M. putorius). We used a juvenile mortality 
rate of 0.7, resulting in a A of 0.5. As with juve¬ 
nile mortality rates, no data on adult mortality 
rates were available for ferrets. Data from other 
mustelids indicated adult mortality ranged from 
15%-25% for Martes pennanti (Kelly 1977), M. 
americana (Strickland et al. 1982), and Mustela 
erminea (Stroganov 1937), although King (1980) 
reported mortality rates of 80%-90% for M. ni¬ 
valis. We varied adult mortality (|x) for ferrets 
from 0.2 to 0.4. 

The probability of ferret populations of size n 
reaching a size where the probability of extinc¬ 
tion is nearly zero was calculated by 
P ~ 1 - (|x/A) n 

where jx and A are the per capita death and birth 
rates, respectively (MacArthur and Wilson 
1967, Richter-Dyn and Goel 1972). Two results 
emerged from this analysis (Table 1). First, as 
jx/A decreases, the number of female ferrets 
needed to colonize an area successfully de¬ 
creases for a given probability. Second, as the 
number of breeding females increases, a higher 
jjl/A can still be tolerated with a high probability 
of successful colonization. These results have a 
practical application in the reintroduction of 
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Table 1. Probability of a female black-footed ferret 
population of size n successfully colonizing and reaching 
carrying capacity. For a birth rate (X) of 0.5 and a death 
rate (|x) varied from 0.2 to 0.4, the resulting (jl/X ranges 
from 0.4 to 0.8. 


Population size (n) Death rate/birth rate ratio (p./X) 
0.9 0.8 0.7 0.6 0.5 0.4 


5 

0.41 

0.67 

0.83 

0.92 0.97 0.99 

10 

0.65 

0.89 

0.97 

0.99 

15 

0.79 

0.96 

0.99 


20 

0.88 

0.99 



30 

0.96 




40 

0.98 




50 

0.99 





black-footed ferrets. For example, if jjl/X is con¬ 
servatively estimated to be 0.9, these data sug¬ 
gest that at least 15 females will be needed for a 
reintroduction that has an 80% chance of suc¬ 
ceeding. As more accurate birth and death rates 
are estimated for the Wyoming ferrets, this type 
of analysis will become even more useful for 
assessing the number of animals needed for 
translocations. 

We also calculated extinction times for an es¬ 
tablished population of ferrets (Table 2). These 
values (T k ) were estimated by the MacArthur 
and Wilson (1967) equations 

T k ~ Tj b where 
X - p. 

Tj = 2 (A/p,)* • (1/iX) 

i- 1 

As the death rate to birth rate ratio (|x/X) in¬ 
creases, the time to extinction (T k ) decreases. 
For a given |x/X, T k increases as the carrying 
capacity (K) increases. Assuming that the carry¬ 
ing capacity of the Wyoming ferret population is 
approximately 40 adults (Forrest et al. 1985) and 
that |x/X equals 0.8, the time to extinction for this 
population as predicted by the T k model would 
be about 1,000 years. Data from Table 1 show 
that this population of 40 (20 females, 20 males) 
with a death rate/birth rate ratio of 0.8 would 
have about a 99% chance of successful coloniza¬ 
tion and conversely less than a 1% probability of 
extinction. Thus, a population of 40 adult ferrets 
should be a reasonable first estimate of MVP size 
as predicted by the T k model. If we assume that 
there is a constant average density of ferrets 
across the study area of 1/50 ha (Forrest et al. 
1985), this MVP estimate would result in a mini¬ 
mum area requirement of 2,000 ha. In actuality, 


Table 2. Times to extinction (T k ) for black-footed fer¬ 
ret populations with a carrying capacity (K) of 40 and 50, 
and birth (X)/ death (|x) rates varied as shown. See text for 
details. 




X 

M- 

Ijl/X 

T k (years) 

k 

40 

0.5 

0.2 

0.4 

1.0 x 10 15 



0.5 

0.3 

0.6 

3.0 x 10 7 



0.5 

0.4 

0.8 

1.0 x 10 3 

K 

50 

0.5 

0.2 

0.4 

8.0 x 10 1S 



0.5 

0.3 

0.6 

5.0 x 10 9 



0.5 

0.4 

0.8 

5.0 x 10 4 


because prairie dog densities vary in this habitat 
(Clark et al. 1985), the 2,000 ha figure is low. 

For endangered species with insular or patchy 
distributions, the theory of island biogeography 
can be a useful tool for estimating MVP size if 
data on birth and death rates are available. The 
T k model is a mathematically simple one that 
assumes limited exponential growth, either 
birth or death rates that vary linearly with den¬ 
sity, no age- specific variance of birth and death 
rates, and a constant environment. Neverthe¬ 
less, the T k model has at least two serious draw¬ 
backs: (1) it ignores inbreeding depression (see 
genetic considerations), which may be impor¬ 
tant in small populations, and (2) violation of the 
model assumptions decreases the time to extinc¬ 
tion (T k ), thus producing what may be overly 
optimistic T k values (Shaffer and Samson 1985). 
For the ferret, there is good evidence that not all 
females breed and that natality varies with age 
(Forrest et al., in manuscript), thereby violating 
two assumptions of the model. The degree to 
which these violations may decrease extinction 
times is unknown. 

Other researchers have criticized all aspects of 
island biogeography theory, the T k model in¬ 
cluded, claiming it remains essentially unsub¬ 
stantiated (Gilbert 1980, Margules et al. 1982). 
Despite these criticisms, several studies have 
shown good agreement between model predic¬ 
tions and field observations. Crowell (1973) re¬ 
ported good accordance between observed and 
predicted times to extinction for mice intro¬ 
duced to islands. Similarly, Smith (1974, 1980) 
and Fritz (1979) predicted times to extinction 
that appeared to explain satisfactorily local dis¬ 
tributions of pikas (Ochotona princeps) and 
spruce grouse (Canachites canadensis). We sug¬ 
gest that the T k model of island biogeography 
remains a useful one that should be used in con¬ 
junction with other methods to estimate MVP 
size and area requirements. 
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Simulation Models 

Computer simulations may be the most useful 
method to estimate MVP size and minimum area 
requirements. Except for theoretical models 
like the T k of island biogeography, computer 
models provide the only method in which a 
probability value for survival can be attached to 
the MVP estimate. In addition, these models are 
not confined to the mathematical assumptions of 
analytical models (e.g., density dependent na¬ 
tality or mortality rates, lack of age structure) and 
provide a flexible mechanism for assessing the 
sensitivity of MVP estimates to changes in cer¬ 
tain population parameters. 

Shaffer (1978) employed the simulation ap¬ 
proach to estimate MVP size and minimum area 
requirements for the grizzly bear (Ursus arctos) 
in Yellowstone National Park. This simulation 
evaluated the effects of both demographic and 
environmental stochasticity and also indicated 
which population parameters were most likely to 
affect changes in survival probability. Watts and 
Conley (1981) used both stochastic and deter¬ 
ministic models to predict survival and extinc¬ 
tion probabilities for a remnant population of 
bighorn sheep in the southwestern United 
States. Their simulations evaluated the effects of 
demographic but not environmental stochasti¬ 
city. Although an estimate of MVP size was not 
an explicit result of their effort, it could have 
been obtained from their simulations. 

The major disadvantage of the simulation ap¬ 
proach is the extensive population data it re¬ 
quires. Minimum data requirements are the 
mean and variance of age-specific and sex- 
specific natality/mortality rates, age structure, 
sex ratio, and the relationship of these variables 
to density (Shaffer 1981). For most species, ob¬ 
taining these data requires an extensive and in¬ 
tensive mark/recapture and/or radio telemetry 
effort over several years. In his grizzly bear sim¬ 
ulation, Shaffer (1978) used the 12-year data 
base of Craighead et al. (1974). Such extensive 
data for any wildlife species are the exception, 
not the rule. However, for some species in which 
data are lacking, it should be feasible to substi¬ 
tute data from a closely related (congeneric) spe¬ 
cies into the simulation. 

The black-footed ferret is a case in point. Al¬ 
though some natality and sex ratio data are avail¬ 
able (Forrest et al., in manuscript), mortality 
data are not. With ferret data on natality and 
steppe polecat (M. eversmanni) data on mortal¬ 


ity, Shaffer (personal communication) simulated 
souk* preliminary estimates of MVP size for the 
ferret. Work is now in progress to refine these 
simulations, which will take into account both 
demographic and environmental stochasticity. 
As the Wyoming ferret study proceeds, more 
age- and sex-specific natality/mortality data may 
become available to incorporate into the simula¬ 
tion. For short-term management needs, how¬ 
ever, a MVP estimate by simulation using data 
from a closely related species should suffice for 
the black-footed ferret. 

Although not directly useful in estimating 
MVP size, simulation models based on bioen¬ 
ergetics are useful in estimating minimum area 
requirements. Two such models exist for the 
ferret. One model estimates gestation, lactation, 
and growth energy requirements for one female 
ferret and her young (Stromberg et al. 1983); the 
other estimates energy requirements based on 
experimental feeding studies of steppe polecats 
and observed activity patterns of the Wyoming 
ferrets (Powell et al. 1985). Both models indicate 
the number of prairie dogs required to support a 
given population of ferrets. Combined with esti¬ 
mates of MVP size, data from these models are 
helpful in estimating minimum area require¬ 
ments for the ferret, assuming that prey (prairie 
dog) densities can be measured and that there is 
a close correspondence between prey density 
and availability. 

Computer simulations, though realistic in 
some aspects, cannot incorporate all the behav¬ 
ioral and ecological adaptations of every species 
(Watts and Conley 1981). They should not be 
used to predict actual population parameters but 
instead should indicate a range of possibilities for 
those parameters (Fowler 1981). In discussing 
the role of computer simulations in wildlife sci¬ 
ence, Romesburg (1981) viewed these models as 
comprehensive tools for integrating knowledge, 
common sense, hunches, and opinions. It is in 
this role that we feel computer simulations can 
be an aid to estimating MVP size and minimum 
area requirements for endangered species. Shaf¬ 
fer’s grizzly bear simulation is presently being 
used by wildlife managers in Montana and Wyo¬ 
ming in this capacity. In addition, an adaptation 
of his simulation is currently being made avail¬ 
able to the US DA Forest Service for use in the 
management of vertebrate species in their forest 
planning process (Shaffer, personal communica¬ 
tion). 
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Genetic Considerations 

Increasing attention is being given to estimat¬ 
ing minimum population numbers from a ge¬ 
netic standpoint (Soule and Wilcox 1980, 
Frankel and Soule 1981, Schonewald-Cox et al. 

1983, Lacava and Hughes 1984, Lehmkuhl 

1984, Salwasser et al. 1984). The critical ques¬ 
tion is what population thresholds are necessary 
to maintain short-term and long-term (evolu¬ 
tionary) fitness. Based on the extensive experi¬ 
ence of animal breeders, Franklin (1980) and 
Soule (1980) determined that the maximum al¬ 
lowable rate of inbreeding necessary to avoid 
short-term effects of inbreeding depression is 
1% per generation. This rule of thumb translates 
to a genetically effective size of 50 for mainte¬ 
nance of short-term fitness. 

The concept of a genetically effective popula¬ 
tion (N e ) is important. Such a population is de¬ 
fined as one in which all individuals mate ran¬ 
domly, the sex ratio approaches unity, variance 
in family size is zero, and generations do not 
overlap. Obviously, few if any vertebrate popu¬ 
lations meet these criteria. Thus, a genetically 
effective size of 50 translates into a larger actual 
(census) number of breeding adults for most spe¬ 
cies. 

Several mathematical expressions have been 
derived for determining the effects of variance in 
progeny, unequal sex ratios, overlapping gener¬ 
ations, and population fluctuations on the genet¬ 
ically effective population size (e.g., Kimura and 
Crow 1963, Emigh and Poliak 1979). Lehmkuhl 
(1984) developed a procedure, based on the 
above formulae, for adjusting the genetically ef¬ 
fective population of 50 to arrive at the census 
number of breeding animals. With Lehmkuhl’s 
procedure, we estimated the number of breed¬ 
ing ferrets (N) necessary to maintain a geneti¬ 
cally effective population (N e ) of 50 (see the Ap¬ 
pendix for calculations). Results of this exercise 
indicated that a MVP of 214 breeding ferrets is 
necessary for maintenance of short-term genetic 
fitness. This MVP estimate is five times the 
known number of adults in the Wyoming popu¬ 
lation (Forrest et al., in manuscript). Minimum 
area requirements for this MVP estimate would 
be about 10,700 ha of white-tailed prairie dog 
habitat. 

A question arises as to what constitutes short¬ 
term versus long-term genetic fitness. This 
question can be satisfactorily answered with the 
following equation 


r=i-(i-_ L)‘ 

2 N e 

where / equals the inbreeding coefficient, N e 
equals the genetically effective population size, 
and t equals generation time. Animal breeders 
have noted a significant reduction in fecundity 
when/approaches0.5-0.6(Soule 1980). Setting 
N e at 50 and/at 0.6 in the above equation results 
in t equal to about 90 generations. If the genera¬ 
tion time for the ferret is 1 year (time from birth 
of a female kit to birth of her first litter), it would 
take about 90 years for an effective population of 
50 ferrets to reach an inbreeding coefficient that 
could cause extinction solely on genetic 
grounds. If the N e for the Wyoming ferret popu¬ 
lation is now substantially less than 50 animals, 
this short-term genetic threshold could be 
reached in considerably fewer than 90 years. 
The amount of time the population has already 
been isolated may be several decades. 

Franklin (1980) suggested that a genetically 
effective population of 500 animals is necessary 
for long-term genetic variation required by the 
evolutionary process. As with short-term fitness, 
a N e of 500 translates into a much larger number 
of breeding adults when dealing with real, as 
opposed to ideal, populations. As an alternative 
to this rule of thumb approach, Frankel and 
Soule (1981) suggested monitoring genetic varia¬ 
tion in target species by electrophoretic tech¬ 
niques that estimate percentage polymorphisms 
and percentage heterozygosity in a population. 
Data from natural populations suggest that rela¬ 
tively heterozygous individuals have greater vi¬ 
ability and fecundity than individuals with lower 
percentages of heterozygosity and polymor¬ 
phism (Soule 1980). 

For many species, particularly endangered 
ones, electrophoretic techniques can be imprac¬ 
tical for several reasons. First, they require po¬ 
tentially injurious tissue sampling. Second, 
small sample sizes of most studies can verify the 
presence but not the absence of rare alleles and 
preclude the use of statistics to test departures 
from Hardy-Weinburg frequencies. Third, 
many higher vertebrates lack polymorphic ge¬ 
netic markers that these techniques can detect. 
Nevertheless, if these techniques become more 
efficient and automated in the future, they could 
be a valuable management tool for working with 
threatened and endangered species. If test re¬ 
sults indicated a reduction in genetic variability, 
then restorative steps, such as increasing N e , 
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could be taken (Frankel and Soule 1981). Base¬ 
line data on genetic variation in the Wyoming 
ferret population are now being gathered (Kil¬ 
patrick et al., 1986). 

Although both empirical data and theory have 
contributed to developing the genetic approach 
to estimating MVP size, this method is not with¬ 
out its shortcomings. Both the “50 and 500 rules” 
which are being advanced are based on simple, 
analytical models of population genetics that do 
not account for either age structure or environ¬ 
mental stochasticity. Previous work indicates 
that the exclusion of age structure can dramati¬ 
cally decrease MVP size (Shaffer 1978) and in so 
doing optimistically deflate the MVP estimates. 

Discussion and Conclusion 

At the outset of this paper, we indicated that 
the five methods for estimating MVP size would 
be evaluated in relation to research impacts on 
target species and the value of accurate and pre¬ 
cise data when time is a critical factor. In the 
black-footed ferret study, we adopted the prag¬ 
matic philosophy of Frankel and Soule (1981) 
that conservationists must often employ rough 
approximations of critical parameters instead of 
waiting for precise data that may never be ob¬ 
tained. Assuming that this philosophy is sound 
for most endangered species studies, in practice 
it leads us to conclude that the experimental and 
simulation methods for estimating MVP size 
have serious limitations. The former requires a 
lengthy study period, and the latter an extensive 
set of population data. Neither time nor exten¬ 
sive population data are goods in great supply for 
most endangered species programs. As previ¬ 
ously pointed out, though, it should be possible 
at times to use less than “perfect” data in the 
simulation approach to produce a meaningful 
MVP estimate. 

Although Shaffer (1981) underscored the im¬ 
portance of defining MVP size with explicit but 
flexible criteria (i.e., time frame and survival 
probability), only two of the five methods exam¬ 
ined (T k and simulation) yield results that in¬ 
clude both of these criteria. Are we then to dis¬ 
miss the other approaches as meaningless 
exercises? We think not. When applicable, both 
the experimental and biogeographical methods 
provide empirical data for estimating MVP size, 
although no survival probability can be attached 
to their estimates. A prudent path for any con¬ 


servation biologist should be to use both empiri¬ 
cal as well as theoretical approaches such as is¬ 
land biogeography and genetic considerations. 

Each of the methods for estimating MVP size 
relies on several assumptions. Which method(s) 
then are most appropriate? Looking at this same 
question from another angle, one might ask to 
which type of perturbation is the target species 
or population most susceptible. For most endan¬ 
gered species, it is unlikely that the latter ques¬ 
tion will be answerable; it is likely that the target 
species will be subject to some combination of 
demographic, environmental, and genetic per¬ 
turbations. Therefore, those methods that take 
into account these perturbations should proba¬ 
bly carry the most weight. Obviously, the avail¬ 
able data will also dictate to some extent which 
methods can be used. 

For the majority of endangered species 
studies, empirical data necessary for estimat¬ 
ing MVP size with the experimental and bio¬ 
geographical methods will not be available. If 
population data on natality and mortality are 
obtainable, then a combination of the simula¬ 
tion and genetic approaches may yield the 
best estimate of MVP size. The simulation 
method can account for both demographic 
and environmental perturbations, and the ge¬ 
netic method can account for demographic 
stochasticity as it relates to genetic drift. On 
the other hand, theoretical methods like the 
T k model of island biogeography only account 
for demographic perturbations. A combina¬ 
tion of the simulation and genetic approaches 
will likely produce a range of MVP estimates. 
Results of the recent Forest Service workshop 
indicated that environmental and demo¬ 
graphic variability may be more important 
than genetic variability in setting the lower 
limit to population viability (M. L. Shaffer, 
personal communication). One possible im¬ 
plication of these results is that the final MVP 
estimate should be biased toward that esti¬ 
mate produced by the simulation method. 
However, management constraints will also 
be a primary factor in establishing MVP size 
somewhere within the range of estimates. 

When no population data are available and 
empirical approaches have also been ruled out, 
then genetic considerations must take priority in 
estimating MVP size. This is the case for the 
ferret. Until more data are obtained for simula¬ 
tions, we must rely on our estimate of about 200 
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animals by the genetics method to serve as the 
MVP estimate for black-footed ferrets. In the 
final analysis, this number may be reduced as 
a result of compromising between estimates 
from the simulation and genetics methods. 

The preceding discussion on genetic ap¬ 
proaches is relevant only to maintenance of 
short-term fitness. For fitness on a long-term 
scale, we must try to establish several ferret 
populations, each of which maintains a short¬ 
term MVP size. Although such goals may now 
appear beyond the vision of what resource 
managers can do to conserve the ferrets in the 
near term, we must never lose sight of the fact 
that large, viable populations of ferrets and 
other endangered species are necessary for 
the recovery and conservation of these species 
from an evolutionary standpoint. 
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Appendix 

Calculation of actual adult census number 
(N) for black-footed ferrets through adjust¬ 
ment of genetically effective population size 
(N e ) (Lehmkuhl 1984). 

Step 1 : Adjustment for variation in progeny number. 

N - (2N e V k + 1)/K 2 

where N e = effective population size of 50, V k variance 
in number of young, K mean number of young per 
female. When lacking data on V k , multiply 1.4 times N t .. 
For the black-footed ferret: 50 x 1.4 70. 


Step 2 : Adjustment for unequal sex ratio. 

N m = (N e + [male:female ratio x N e ])/4 
where N m = number of males. 

Nf = N m X female:male ratio 

N = N m + N f 

N/50 = increase in N e 

N = increase X 70 (from Step 1) 

For the black-footed ferret: (sex ratio data from 1984 
census—Forrest et al., in manuscript) 

N m - (50 + 18/25(50))/4 - 21.5 
Nf - 21.5 x 25/18 - 29.9 
N = 29.9 + 21.5 = 51.4 
51.4/50 - 1.02 
1.02 X 70 - 71.4 

Step 3 : Adjustment for overlapping generation. 

Double the census number from Step 2. 

For the black-footed ferret: 2x 71.4 =142.8 
Step 4 : Adjustment for population fluctuations. 

Ratio between empirical high and low population cen¬ 
suses indicates factor with which to increase Step 3 result 
for MVP estimate. 

For black-footed ferret: 43/28 1.5 1.5 X 142.8 = 214 

MVP (1983 census = 28 adults, 1984 = 43, Forrest et 
al., in manuscript) 
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